System and method for tracking body joints

ABSTRACT

Body joint tracking is applied in various industries and medical field. In body joint tracking, marker less devices plays an important role. However, the marker less devices are facing some challenges in providing optimal tracking due to occlusion, ambiguity, lighting conditions, dynamic objects etc. System and method of the present disclosure provides an optimized body joint tracking. Here, motion data pertaining to a first set of motion frames from a motion sensor are received. Further, the motion data are processed to obtain a plurality of 3 dimensional cylindrical models. Here, every cylindrical model among the plurality of 3 dimensional cylindrical model represents a body segment. The coefficients associated with the plurality of 3 dimensional cylindrical models are initialized to obtain a set of initialized cylindrical models. A set of dynamic coefficients associated with the initialized cylindrical models are utilized to track joint motion trajectories of a set of subsequent frames.

PRIORITY CLAIM

This U.S. patent application claims priority under 35 U.S.C. § 119 to: India Application No. 201721024066, filed on Jul. 7, 2017. The entire contents of the aforementioned application are incorporated herein by reference.

TECHNICAL FIELD

The embodiments herein generally relates, in general, to motion tracking and, in particular, to a system and method for tracking body joints.

BACKGROUND

Tracking body joints is key research field and applied in a variety of applications including computer animation, video games, medical therapy, surveillance, human machine interaction and athlete performance analysis. Moreover, the joint motion analysis is an important aspect in health monitoring of patients suffering from neurological disorders, post-stroke patients, and an elderly subjects. In real-time body joint tracking, accuracy in tracking is a challenging task. The challenge is due to change in appearance due to non-rigid motion of a subject undergoing test, clothing, view point and lighting.

Conventional real-time body joint tracking methods are marker-based and requires a subject undergoing test to wear obtrusive devices. Moreover, the conventional marker based joint tracking methods are complex, difficult to maintain, extremely expensive and may not be viable for a prolonged rehabilitation therapy. Hence marker less joint tracking systems are used for tracking body joints in real-time. However, the conventional marker less real-time body joint tracking systems are facing challenges to obtain optimized body joint tracking due to occlusion, ambiguity, lighting conditions and dynamic objects and the like.

SUMMARY

Embodiments of the present disclosure present technological improvements as solutions to one or more of the above-mentioned technical problems recognized by the inventors in conventional systems. For example, in one embodiment, a method for tracking body joints is provided. The method includes receiving, motion data pertaining to a first set of motion frames from a motion sensor device, wherein the motion data comprises an initial joint motion data and a plurality of depth frames, wherein the initial joint motion data comprises locations in 3 dimensional world coordinates of a plurality of body joints associated with the plurality of depth frames and each depth frame among the plurality of depth frames comprises a set of 2 dimensional image coordinates and a depth value, by one or more hardware processors. Further, the method includes, initializing, a plurality of cylindrical models for each motion frame among the first set of motion frames by processing the initial joint motion data and the plurality of depth frames, wherein the plurality of cylindrical models represents a plurality of body segments, connecting adjacent body joints among the plurality of body joints, by the one or more hardware processors. Furthermore, the method includes tracking, by the one or more hardware processors, the plurality of initialized cylindrical models to obtain one or more optimized motion trajectories of the plurality of body joints in a plurality of motion frames received in succession to the first set of motion frames based on one or more direction angles and one or more base coordinates associated with each of the plurality of cylindrical models, wherein the tracking is performed by utilizing particle filter mechanism, by the one or more hardware processors.

In another aspect, a system for tracking body joints is provided. The system includes one or more memories comprising programmed instructions and repository for storing the one or more motion data, a database; and one or more hardware processors operatively coupled to the one or more memories, wherein the one or more hardware processors are capable of executing the programmed instructions stored in the one or more memories, a motion sensor device and a joint tracking unit, wherein the joint tracking unit is configured to receive, motion data pertaining to a first set of motion frames from a motion sensor device, wherein the motion data comprises an initial joint motion data and a plurality of depth frames, wherein the initial joint motion data comprises locations in 3 dimensional world coordinates of a plurality of body joints associated with the plurality of depth frames and each depth frame among the plurality of depth frames comprises a set of 2 dimensional image coordinates and a depth value. Further, the joint tracking unit is configured to initialize, a plurality of cylindrical models for each motion frame among the first set of motion frames by processing the initial joint motion data and the plurality of depth frames, wherein the plurality of cylindrical models represents a plurality of body segments, connecting adjacent body joints among the plurality of body joints. Furthermore, the joint tracking unit is configured to track, the plurality of initialized cylindrical models to obtain one or more optimized motion trajectories of the plurality of body joints in a plurality of motion frames received in succession to the first set of motion frames based on one or more direction angles and one or more base coordinates associated with each of the plurality of cylindrical models, wherein the tracking is performed by utilizing particle filter mechanism.

In yet another aspect, a computer program product comprising a non-transitory computer-readable medium having embodied therein a computer program for system and method for tracking body joints, is provided. The computer readable program, when executed on a computing device, causes the computing device to receive, motion data pertaining to a first set of motion frames from a motion sensor device, wherein the motion data comprises an initial joint motion data and a plurality of depth frames, wherein the initial joint motion data comprises locations in 3 dimensional world coordinates of a plurality of body joints associated with the plurality of depth frames and each depth frame among the plurality of depth frames comprises a set of 2 dimensional image coordinates and a depth value. Further computer readable program, when executed on a computing device, causes the computing device to initialize, a plurality of cylindrical models for each motion frame among the first set of motion frames by processing the initial joint motion data and the plurality of depth frames, wherein the plurality of cylindrical models represents a plurality of body segments, connecting adjacent body joints among the plurality of body joints. Furthermore computer readable program, when executed on a computing device, causes the computing device to track, the plurality of initialized cylindrical models to obtain one or more optimized motion trajectories of the plurality of body joints in a plurality of motion frames received in succession to the first set of motion frames based on one or more direction angles and one or more base coordinates associated with each of the plurality of cylindrical models, wherein the tracking is performed by utilizing particle filter mechanism.

It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the invention, as claimed.

BRIEF DESCRIPTION OF THE DRAWINGS

The accompanying drawings, which are incorporated in and constitute a part of this disclosure, illustrate exemplary embodiments and, together with the description, serve to explain the disclosed principles:

FIG. 1 illustrates a network environment implementing a system and method for tracking body joints, according to some embodiments of the present disclosure;

FIG. 2 illustrates a block diagram of a system and method for tracking body joints, according to some embodiments of the present disclosure;

FIG. 3 depicts an architecture diagram for a system and method for tracking body joints according to some embodiments of the present disclosure;

FIG. 4 illustrates a detailed flow diagram for tracking body joints, according to some embodiments of the present disclosure.

FIG. 5A illustrates an example 3 dimensional temporal trajectory of wrist joint during shoulder abduction, according to some embodiments of the present disclosure;

FIG. 5B illustrates an example 3 dimensional temporal trajectory of elbow joint during shoulder flexion, according to some embodiments of the present disclosure;

and

FIG. 5C illustrates an example 3 dimensional temporal trajectory of wrist joint during elbow flexion, according to some embodiments of the present disclosure.

It should be appreciated by those skilled in the art that any block diagrams herein represent conceptual views of illustrative systems and devices embodying the principles of the present subject matter. Similarly, it will be appreciated that any flow charts, flow diagrams, and the like represent various processes which may be substantially represented in computer readable medium and so executed by a computer or processor, whether or not such computer or processor is explicitly shown.

DETAILED DESCRIPTION

Exemplary embodiments are described with reference to the accompanying drawings. In the figures, the left-most digit(s) of a reference number identifies the figure in which the reference number first appears. Wherever convenient, the same reference numbers are used throughout the drawings to refer to the same or like parts. While examples and features of disclosed principles are described herein, modifications, adaptations, and other implementations are possible without departing from the spirit and scope of the disclosed embodiments. It is intended that the following detailed description be considered as exemplary only, with the true scope and spirit being indicated by the following claims.

In the field of motion tracking, tracking body joints plays a vital role. Since the conventional methods are performing tracking of body joints by utilizing marker based methods and the marker based methods needs some obtrusive devices to be worn by the subjects undergoing tests, marker less body joints tracking methods are prominent. The marker less motion body joint tracking methods includes a motion sensor to capture a plurality of motion data. For example, Kinect® V1 is marker less motion sensing device mainly used for tracking body joints. However, the motion sensor based body joint tracking methods are prone to occlusion, IR interference, lighting conditions and distance from the motion sensor.

The present subject matter overcomes the limitations of the conventional joint tracking methods by modeling a set of 3 dimensional point clouds associated with a body (for example, human body and animal body). Here, a set of 3 dimensional cylindrical models are initialized by utilizing the set of 3 dimensional point cloud data. Further, the set of initialized 3 dimensional cylindrical models are tracked based on a set of dynamic parameters associated with the set of initialized cylindrical models to obtain an optimized joint motion trajectories associated with the body. Here, the optimized joint locations are derived directly in 3 dimensional coordinate space to reduce adverse effects of noise associated with a depth data from depth sensor of a motion sensor device. An implementation of the system and method for tracking body joints is described further in detail with reference to FIGS. 1 through 5C.

Referring now to the drawings, and more particularly to FIGS. 1 through 5C, where similar reference characters denote corresponding features consistently throughout the figures, there are shown preferred embodiments and these embodiments are described in the context of the following exemplary system and/or method.

FIG. 1 illustrates a network environment 100 implementing a system 102 for tracking body joints, according to an example embodiment of the present subject matter. The system for tracking body joints 102, hereinafter referred to as the system 102, is configured for tracking body joints using motion data received from a motion sensor device 120. In an embodiment, motion sensor device 120 can be a Kinect® V1 device for monitoring subject's movements. The system 102 may be embodied in a computing device, for instance a computing device 104.

Although the present disclosure is explained considering that the system 102 is implemented on a server, it may be understood that the system 102 may also be implemented in a variety of computing systems, such as a laptop computer, a desktop computer, a notebook, a workstation, a cloud-based computing environment and the like. In one implementation, the system 102 may be implemented in a cloud-based environment. It will be understood that the system 102 may be accessed by multiple users through one or more user devices 106-1, 106-2 . . . 106-N, collectively referred to as user devices 106 hereinafter, or applications residing on the user devices 106. Examples of the user devices 106 may include, but are not limited to, a portable computer, a personal digital assistant, a handheld device, a Smartphone, a Tablet Computer, a workstation and the like. The user devices 106 are communicatively coupled to the system 102 through a network 108.

In an embodiment, the network 108 may be a wireless or a wired network, or a combination thereof. In an example, the network 108 can be implemented as a computer network, as one of the different types of networks, such as virtual private network (VPN), intranet, local area network (LAN), wide area network (WAN), the internet, and such. The network 106 may either be a dedicated network or a shared network, which represents an association of the different types of networks that use a variety of protocols, for example, Hypertext Transfer Protocol (HTTP), Transmission Control Protocol/Internet Protocol (TCP/IP), and Wireless Application Protocol (WAP), to communicate with each other. Further, the network 108 may include a variety of network devices, including routers, bridges, servers, computing devices, storage devices. The network devices within the network 108 may interact with the system 102 through communication links.

As discussed above, the system 102 may be implemented in a computing device 104, such as a hand-held device, a laptop or other portable computer, a tablet computer, a mobile phone, a PDA, a smartphone, and a desktop computer. The system 102 may also be implemented in a workstation, a mainframe computer, a server, and a network server. In an embodiment, the system 102 may be coupled to a data repository, for example, a repository 112. The repository 112 may store data processed, received, and generated by the system 102. In an alternate embodiment, the system 102 may include the data repository 112. The components and functionalities of the system 102 are described further in detail with reference to FIG. 2.

FIG. 2 illustrates a block diagram of a body joint tracking system 200 for tracking body joints, according to some embodiments of the present disclosure. The body joint tracking system 200 (hereinafter referred to as system 200) may be an example of the system 102 (FIG. 1). In an example embodiment, the system 200 may be embodied in, or is in direct communication with the system, for example the system 102 (FIG. 1). The system 200 includes or is otherwise in communication with one or more hardware processors such as a processor 202, at least one memory such as a memory 204, an I/O interface 206 and a joint tracking unit 250. In an embodiment, the joint tracking unit 250 can be implemented as a standalone unit in the system 200 comprising an point cloud segmentation module (not shown in FIG. 2), a cylinder model fitting module (not shown in FIG. 2), a state propagation module (not shown in FIG. 2) and a likelihood estimation module (not shown in FIG. 2). In another embodiment, the joint tracking unit 250 can be implemented as a module in the memory 204 comprising the point cloud segmentation module (not shown in FIG. 2), the cylinder model fitting module (not shown in FIG. 2), the state propagation module (not shown in FIG. 2) and the likelihood estimation module (not shown in FIG. 2). The processor 202, memory 204, and the I/O interface 206 may be coupled by a system bus such as a system bus 208 or a similar mechanism.

The I/O interface 206 may include a variety of software and hardware interfaces, for example, a web interface, a graphical user interface, and the like. The interfaces 206 may include a variety of software and hardware interfaces, for example, interfaces for peripheral device(s), such as a keyboard, a mouse, an external memory, a camera device, and a printer. Further, the interfaces 206 may enable the system 102 to communicate with other devices, such as web servers and external databases. The interfaces 206 can facilitate multiple communications within a wide variety of networks and protocol types, including wired networks, for example, local area network (LAN), cable, etc., and wireless networks, such as Wireless LAN (WLAN), cellular, or satellite. For the purpose, the interfaces 206 may include one or more ports for connecting a number of computing systems with one another or to another server computer. The I/O interface 206 may include one or more ports for connecting a number of devices to one another or to another server.

The hardware processor 202 may be implemented as one or more microprocessors, microcomputers, microcontrollers, digital signal processors, central processing units, state machines, logic circuitries, and/or any devices that manipulate signals based on operational instructions. Among other capabilities, the hardware processor 202 is configured to fetch and execute computer-readable instructions stored in the memory 204.

The memory 204 may include any computer-readable medium known in the art including, for example, volatile memory, such as static random access memory (SRAM) and dynamic random access memory (DRAM), and/or non-volatile memory, such as read only memory (ROM), erasable programmable ROM, flash memories, hard disks, optical disks, and magnetic tapes. In an embodiment, the memory 204 includes a plurality of modules 220 and a repository 240 for storing data processed, received, and generated by one or more of the modules 220 and the joint tracking unit 250. The modules 220 may include routines, programs, objects, components, data structures, and so on, which perform particular tasks or implement particular abstract data types.

The memory 204 also includes module(s) 220 and a data repository 240. The module(s) 220 include programs or coded instructions that supplement applications or functions performed by the body joint tracking system 200. The modules 220, amongst other things, can include routines, programs, objects, components, and data structures, which perform particular tasks or implement particular abstract data types. The modules 220 may also be used as, signal processor(s), state machine(s), logic circuitries, and/or any other device or component that manipulates signals based on operational instructions. Further, the modules 220 can be used by hardware, by computer-readable instructions executed by a processing unit, or by a combination thereof. The modules 220 can include various sub-modules (not shown). The module 220 may include computer-readable instructions that supplement applications or functions performed by the signal analysis system 200.

The data repository 240 may include received motion tracking data 242, 3D world coordinate data 244, a plurality of body joints 246 and other data 248. Further, the other data 248 amongst other things, may serve as a repository for storing data that is processed, received, or generated as a result of the execution of one or more modules in the module(s) 220 and the modules associated with the joint tracking unit 250.

Although the data repository 240 is shown internal to the body joint tracking system 200, it will be noted that, in alternate embodiments, the data repository 240 can also be implemented external to the body joint tracking system 200, where the data repository 240 may be stored within a database (not shown in FIG. 2) communicatively coupled to the body joint tracking system 200. The data contained within such external database may be periodically updated. For example, new data may be added into the database (not shown in FIG. 2) and/or existing data may be modified and/or non-useful data may be deleted from the database (not shown in FIG. 2). In one example, the data may be stored in an external system, such as a Lightweight Directory Access Protocol (LDAP) directory and a Relational Database Management System (RDBMS). In another embodiment, the data stored in the data repository 240 may be distributed between the body joint tracking system 200 and the external database.

The joint tracking unit 250 of the body joint tracking system 200 can be configured to receive, a plurality of motion data pertaining to a first set of motion frames from the motion sensor device 120. The motion data comprises an initial joint motion data and a plurality of depth frames, wherein the initial joint motion data comprises locations in 3 dimensional world coordinates represented as {right arrow over (C)}_(t) ^(j)=[X, Y, Z]^(T) of a plurality of body joints associated with the plurality of depth frames and each depth frame among the plurality of depth frames comprises a set of 2D image coordinates and a depth value. The depth value can be a distance in millimeter between the subject undergoing exercise and the motion sensor device 120. A line joining two adjacent body joints i and j can be {right arrow over (b)}_(t) ^((i,j)). The line {right arrow over (b)}_(t) ^((i,j)) can be represented as a body segment B^((i,j)) connecting two adjacent body joints i and j. Further, the line {right arrow over (b)}_(t) ^((i,j)) is associated with a direction information. In an embodiment, the Kinect® V1 device is utilized for capturing the plurality of motion data including the initial joint motion data and the plurality of depth frames at 30 fps, at a distance of 1.8 m to 2.5 m from the subject undergoing exercise.

In an embodiment, a plurality of active joint Range Of Motion (ROM) data pertaining to upper body portion (for example, shoulder abduction, shoulder flexion and extension, Elbow flexion and extension) of the subject undergoing exercise is received from the motion sensor device 120. In the context of present disclosure, ROM refers to any kind of movement of a body during dynamic postures involving joints in the body such as ROM exercises.

Further, the joint tracking unit 250 of the body joint tracking system 200 can be further configured to initialize a plurality of cylindrical models for each motion frame among the first set of motion frames by processing the initial joint motion data and the plurality of depth frames. Here, the plurality of cylindrical models represents a plurality of body segments, connecting adjacent body joints among the plurality of body joints. Here, a plurality of 3 dimensional point clouds corresponding to the plurality of depth frames are obtained by mapping every 2D pixel coordinates (x, y) to the 3 dimensional world coordinates (X,Y,Z) based on a perspective projection as given below:

${X = {\frac{\left( {x - c_{x}} \right)}{f_{x}}d}},{Y = {\frac{\left( {y - c_{y}} \right)}{f_{y}}d}},$

Z=d, where ƒ_(x) and ƒ_(y) refers to focal lengths of a depth sensor associated with the motion sensor device 120 along x axis and y axis respectively, d is the depth value and (c_(x),c_(y)) are a coordinates of an image center. Further, the plurality of 3 dimensional point clouds P_(t) are segmented based on a segmentation threshold to obtain a first set of segmented 3 dimensional point clouds S_(t) representing a body segment B^((i,j)) using the equation,

S _(t) {{right arrow over (P)}∈P _(t):∥({right arrow over (P)}−{right arrow over (b)} _(t) ^((i,j)))∥<r}  →(1)

Here, r is the segmentation threshold obtained based on an Euclidean distance computed for 3 dimensional world coordinates corresponding to the plurality of body joints of the initial joint motion data. Here, the first set of segmented 3 dimensional point clouds S_(t) corresponds to the plurality of body segments B^((i,j)) for every joint pair (i, j). Further, the first set of segmented 3 dimensional point clouds S_(t) corresponding to the plurality of body segments B^((i,j)) are smoothened and reconstructed prior to model fitting to reduce the model outliers. The smoothened first set of segmented 3 dimensional point clouds S′_(t) are reconstructed and smoothened by applying at least one surface reconstruction mechanism. In an embodiment, Moving Least Squares (MLS) can be utilized for the surface reconstruction and smoothening. The surface reconstruction and smoothening on the first set of segmented 3 dimensional point clouds S_(t) provides a new set of segmented point cloud S′_(t) corresponding to the plurality of body segments B^((i,j)). Further, the plurality of cylindrical models are fitted to the new set of segmented point cloud S′_(t) corresponding to the plurality of body segments B^((i,j)). In an embodiment, M_(t) denotes a best fitting cylindrical model at time t with maximum inlier count. Each cylindrical model among the plurality of cylindrical models includes a set of cylinder coefficients represented by M^((i,j))={{right arrow over (b)}, {right arrow over (d)}, R, L}, where, R denotes radius, L denotes length, {right arrow over (d)} denotes direction angles and {right arrow over (b)} denotes base coordinates. The set of cylindrical coefficients M^((i,j)) are obtained from the plurality of cylindrical models M_(t) at different time instances t=1, 2, 3 . . . n, associated with each motion frame among the first set of motion frames. The radius R of the plurality of cylindrical models can be initialized to a mean radius value obtained from the cylinder radius in M_(t) corresponding to the plurality of motion frames. The direction angles {right arrow over (d)} of the plurality of cylindrical models is initialized to a mean direction value {right arrow over (d)}=[d_(x), d_(y), d_(z)]^(T) obtained from the cylinder directions in M_(t) corresponding to from the plurality of motion frames. The base {right arrow over (b)} of the plurality of cylindrical models can be initiated to an average base value {right arrow over (b)}=[b_(x), b_(y), b_(z)]^(T), obtained by projecting a plurality of joint center coordinates {right arrow over (C)}_(t) ^(i) in the initial joint motion data to an axis of a corresponding cylindrical model M_(t) from the plurality of motion frames. The length L of the plurality of cylindrical models can be initialized to a mean length value represented as ∥{right arrow over (b)}_(t) ^((i,j))∥(t=1, 2 . . . n), obtained from a distance between adjacent body joints in the initial joint motion data for plurality of motion frames.

Further the joint tracking unit 250 of the body joint tracking system 200 can be configured to track the plurality of initialized cylindrical models to recover motion trajectories Tr_(j)={({right arrow over (C)}_(t) ^(j), t)|t=1, 2 . . . T} of the plurality of body joints in a plurality of motion frames received in succession to the first set of motion frames based on direction angles {right arrow over (d)} and base {right arrow over (b)} coordinates associated with each of the plurality of cylindrical models. Here, the tracking is performed by utilizing a particle filter mechanism. The method of tracking by utilizing the particle filter mechanism is explained as follows: Here, a plurality of initialized cylindrical models and a second set of segmented point cloud data are received. The second set of segmented point cloud data corresponds to the plurality of motion frames received in succession to the first set of motion frames from the motion sensor device 120. The plurality of initialized cylindrical models are represented as a plurality of particles p^((k)) where k=1, 2 . . . N. The particle filter mechanism approximates a state conditional density of the initialized cylindrical model M^((i,j)) by utilizing a set of N particles p^((k)), where k=1, 2 . . . N. Every particle p^((k)) is associated with a state X^((k)) and a weight π^((k)). The state and the weight values can be updated iteratively. The state X^((k)) of a particle is represented by the direction angles {right arrow over (d)} and base {right arrow over (b)} coordinates, and a weight. At every time instance t, the plurality of particles p^((k)) propagates to a new state. The propagation is determined by a state space model based on a state propagation rule. The state propagation rule is obtained from the plurality of initialized cylindrical models. For example, At time instance t, the particle p^((k)), has a state X_(t) ^((k))=[{right arrow over (b)}_(t), {right arrow over (θ)}_(t)]^(T) where {right arrow over (b)}_(t)=[b_(x), b_(y), b_(z)]^(T) represents the cylinder base coordinates and {right arrow over (θ)}_(t)=[θ_(x), θ_(y), θ_(z)]^(T) represents the direction angles with x, y and z axes. The direction of the cylindrical model represented by particle p^((k)) among the plurality of cylindrical models can be obtained as {right arrow over (d)}_(t)=[d_(x), d_(y), d_(z)]^(T)=[cos θ_(x), cos θ_(y), cos θ_(z)]^(T). The state X_(t) ^((k)) of the particle p^((k)), at a time instance t is updated according to the state propagation rule referred to as equation 2 and equation 3.

{right arrow over (θ)}_(t)={right arrow over (θ)}_(t-1)+({right arrow over (θ)}_(t-1)−{right arrow over (θ)}_(t-2))+N(θ,σ_(θ) ²)  →(2)

{right arrow over (b)} _(t) ={right arrow over (b)} _(t-1) +N(0,σ_(b) ²)  →(3)

Where, σ_(θ) ² and σ_(b) ² are variance of Gaussian noise with zero mean for {right arrow over (θ)} and {right arrow over (b)} respectively. The weight π^((k)) of a particle P^((k)) among the plurality of particles represents a likelihood of observing the new segmented point cloud S′_(t) for a propagated state. The updating of weight π^((k)) is performed by an observation model based on a weight update rule. The weight update rule is obtained from the second set of segmented point cloud data and the plurality of initialized cylindrical models. For example, the weight π^((k)) of the particle P^((k)) can be updated by utilizing the weight update rule referred to as equation 4.

$\begin{matrix} {\pi_{t}^{(k)} = \left. e^{- \frac{D^{2}}{2\; \sigma^{2}}}\rightarrow \right.} & (4) \end{matrix}$

Where σ² is variance of Gaussian likelihood function. D is the sum of all distances d of the cloud points {right arrow over (P)}∈S′_(t) from the cylindrical model represented by a particle p^((k)). The distance d can be calculated by utilizing equation 5.

d=ad _(n)+(1−α)(d _(e) −R)  →(5)

Where, 0<α<1 and d_(n), represents an angular distance between a surface normal {circumflex over (n)}_(p) of the point cloud {right arrow over (P)} and a unit vector {circumflex over (n)}, normal to a cylindrical surface of the cylindrical model and d_(e) is the Euclidean distance between {right arrow over (P)} and the axis of the cylindrical model {right arrow over (d)}_(t). If the distance of new segmented point clouds S′_(t) from the axis of the cylindrical model {right arrow over (d)}_(t) of the particle P^((k)) is lower the probability of correctness of state prediction of the particle is greater. Further, the cylindrical coefficients of the optimal cylindrical model M′^((i,j)). are obtained from a weighted average of the set of cylindrical coefficients—of N number of particles P^((k)). Further, the optimized temporal trajectories Tr_(i) and Tr_(j) corresponding to the joints i and j of the optimal cylindrical model are recovered from a set of end points of the optimal cylindrical model M′^((i,j)).

FIG. 3 depicts an architecture diagram for a system and method for body joint tracking according to some embodiments of the present disclosure. Here, the first set of point cloud data corresponding to the initial set of motion frames and the initial joint motion data are given to the point cloud segmentation module 302. The point cloud segmentation module 302 segments the point cloud data corresponding to the initial set of motion frames into a first set of segmented point clouds. Here the first set of segmented point clouds represents the plurality of body segments. The segmentation is performed by utilizing the initial joint motion data. Further, the plurality of cylindrical models are fitted into the first set of segmented point clouds by the cylinder model fitting module 304. Here, the set of coefficients associated with the plurality of cylindrical models are initialized to obtain the plurality of initialized cylindrical models. Every cylinder among the plurality of initialized cylindrical models represents the corresponding body segment. Here the body segment is connecting two adjacent body joints. Further, a second set of segmented point cloud data corresponding to the plurality of motion frames received in succession to the first set of motion frames from the motion sensor device is received. The second set of point cloud data is segmented using the initial joint motion data by the point cloud segmentation module 310. Further, everybody segment corresponding to the plurality of initialized cylindrical models are represented as the particle. Every particle is associated with the corresponding state and the corresponding weight. The state includes base of the initialized cylindrical model and the direction angles of the initialized cylindrical model. The states are propagated based on the state propagation rule by the state propagation module 306. The state propagation module 306 utilizes the state space model and the state propagation rule is based on the plurality of initialized cylindrical models. The weight associated with the particle is updated by the likelihood estimation module 308 based on the weight update rule. The weight update rule is based on the second set of segmented point cloud data and the initialized cylindrical model. The likelihood estimation module 308 utilizes the observation module. Further, a cylindrical model among the plurality of initialized cylindrical model having maximum weight is selected as an optimal cylindrical model. Here, the cylindrical coefficients of the optimal cylindrical model are obtained from a weighted average of the set of cylindrical coefficients of N number of particles associated with the plurality of initialized cylindrical models. Further, the optimized temporal trajectories Tr_(i) and Tr_(j) corresponding to the joints i and j of the optimal cylindrical model are recovered from a set of end points of the optimal cylindrical model.

FIG. 4 illustrates a flow diagram of a method 400 for body joint tracking, according to some embodiments of the present disclosure. The method 400 may be described in the general context of computer executable instructions. Generally, computer executable instructions can include routines, programs, objects, components, data structures, procedures, modules, functions, etc., that perform particular functions or implement particular abstract data types. The method 400 may also be practiced in a distributed computing environment where functions are performed by remote processing devices that are linked through a communication network. The order in which the method 400 is described is not intended to be construed as a limitation, and any number of the described method blocks can be combined in any order to implement the method 400, or an alternative method. Furthermore, the method 400 can be implemented in any suitable hardware, software, firmware, or combination thereof.

At 402, the system 200 receives, by the one or more hardware processors, motion data pertaining to a first set of motion frames from the motion sensor device. Here, the motion data comprises the initial joint motion data and a plurality of depth frames. The initial joint motion data comprises 3 dimensional world coordinates of the plurality of body joints associated with the plurality of depth frame and each depth frame among the plurality of depth frames comprises the set of 2 dimensional image coordinates and the depth value. At 404, the system 200 initializes a plurality of cylindrical models, by the one or more hardware processors, for each motion frame among the first set of motion frames by processing the initial joint motion data and the plurality of depth frames. Here, the plurality of cylindrical models represents the plurality of body segments, connecting adjacent body joints among the plurality of body joints. Further, the plurality of 3 dimensional point clouds corresponding to the plurality of depth frames are obtained by mapping every 2 dimensional image coordinates to the 3 dimensional world coordinates. The plurality of 3 dimensional point clouds are segmented based on a segmentation threshold to obtain a first set of segmented 3 dimensional point clouds. Here, the first set of segmented 3 dimensional point clouds corresponds to the plurality of body segments. Additionally, the segmentation threshold is obtained based on a Euclidean distance computed for 3 dimensional world coordinates corresponding to the plurality of body joints of the initial joint motion data. Further, the first set of segmented 3 dimensional point clouds corresponding to the plurality of body segments are smoothened prior to model fitting to reduce one or more model outliers. Further, the plurality of cylindrical models are fitted to the set of smoothened point clouds corresponding to the plurality of body segments associated with each motion frame among the first set of motion frames. The set of cylinder coefficients of each cylindrical model among the plurality of cylindrical models includes radius, length, one or more direction angles and one or more base coefficients. The radius of the plurality of cylindrical models are initialized to a mean radius value obtained from the cylinder radius corresponding to the plurality of motion frames. The one or more direction angles of the plurality of cylindrical models are initialized to a mean direction value obtained from a cylinder directions corresponding to the plurality of motion frames. The one or more base coordinates of the plurality of cylindrical models are initialized to an average base value obtained by projecting a plurality of joint center coordinates in the initial joint motion data to an axis of a corresponding cylindrical model from the plurality of motion frames. The length of the plurality of cylindrical models are initialized to a mean length value obtained from a distance between adjacent body joints in the initial joint motion data for the plurality of motion frames. At 406, the system 200 tracks, by the one or more hardware processors, the plurality of initialized cylindrical models to obtain one or more optimized motion trajectories of the plurality of body joints in a plurality of motion frames received in succession to the first set of motion frames based on the one or more direction angles and the one or more base coordinates associated with each of the plurality of cylindrical models. Here, the tracking is performed by utilizing the particle filter mechanism, the particle filter mechanism receives the plurality of initialized cylindrical models and the second set of segmented point cloud data corresponding to the plurality of motion frames received in succession to the first set of motion frames from the motion sensor device. Here, the initialized plurality of cylindrical models are represented as a plurality of particles. Each particle among the plurality of particles comprises a state represented by the one or more direction and the one or more base coordinates and a weight. The state of the plurality of particles are propagated to a new state and the state is updated based on the state propagation rule by utilizing the state space model, wherein the state propagation rule is based on the plurality of initialized cylindrical models. The weight associated with the plurality of particles are updated based on the likelihood function and the weight update rule by utilizing an observation model, wherein the weight update rule is based on the second set of segmented point cloud data and the plurality of initialized cylindrical models. Further, the plurality of initialized cylindrical models are tacked based on the state and the weight of the plurality of particles associated with the plurality of initialized cylindrical models to obtain the optimized motion trajectories of the plurality of body joints.

In an embodiment, the system 200 is experimented as follows: The temporal trajectory of a joint among the plurality of joints moving with respect to a stationary pivot joint (e.g trajectory of the elbow about the shoulder joint during shoulder abduction) is analyzed for motion smoothness and variation of radius. FIG. 5A to 5C illustrates examples of 3 dimensional trajectory plots comparing smoother trajectory of motion and lesser variation in radius of the trajectories obtained after applying the present disclosure, with Kinect® V1 skeletal tracking. Particularly, FIG. 5A illustrates an example 3 dimensional temporal trajectory of wrist joint during shoulder abduction, according to some embodiments of the present disclosure. Here, a plot with black circles illustrates the temporal trajectory of the wrist obtained through the present disclosure and a plot with white squares illustrates the temporal trajectory of the wrist obtained through Kinect® V1. FIG. 5B illustrates an example 3 dimensional temporal trajectory of elbow joint during shoulder flexion, according to some embodiments of the present disclosure. Here, the plot with black circles illustrates the temporal trajectory of elbow joint during shoulder flexion obtained through the present disclosure and the plot with white squares illustrates the temporal trajectory of elbow joint during shoulder flexion obtained through Kinect® V1. FIG. 5C illustrates an example 3 dimensional temporal trajectory of wrist joint during elbow flexion, according to some embodiments of the present disclosure. Here, the plot with black circles illustrates the temporal temporal trajectory of wrist joint during elbow flexion obtained through the present disclosure and the plot with white squares illustrates the temporal temporal trajectory of wrist joint during elbow flexion obtained through Kinect®

In an embodiment, the following performance metrics are computed: (i) root-mean-squared error (RMSE), (ii) coefficient of variation (CV) of the radius of curvature, and (iii) a measure for motion smoothness. The RMSE is computed with respect to a reference radius of curvature (equal to the body segment length) estimated during initialization of the plurality of cylinders. The coefficient of variation (CV) is a measure of trajectory spread that describes the ratio of variance to mean of the trajectory radius. The smoothness of a function ƒ(x) approximating the trajectory points over an interval [a;b] is measured by function g(ƒ(x)) referred to as equation 6, mainly used for smoothing spline interpolation. A lower value of g(ƒ(x)) indicates greater smoothness of ƒ(x).

g(ƒ(x))=∫_(a) ^(b)(ƒ″(x))² dx  →(6)

TABLE I ROM Trajectory Coefficient RMSE Smoothness Exercise of joint Method of variation (in mm) coefficient Shoulder Elbow Kinect ® V1 0.128870 0.051360 0.063295 Abduction- Present 0.091506 0.042277 0.043249 Adduction disclosure Wrist Kinect ® V1 0.084366 0.130863 0.184191 Present 0.053918 0.044017 0.077605 disclosure Shoulder Elbow Kinect ® V1 0.101716 0.039362 0.039370 Flexion- Present 0.070454 0.028332 0.028636 Extension disclosure Wrist Kinect ® V1 0.060563 0.128109 0.092172 Present 0.042490 0.028614 0.043032 disclosure Elbow Wrist Kinect ® V1 0.069170 0.086471 0.017815 Flexion- Present 0.038571 0.023394 0.016184 Extension disclosure Percentage of Improvement 33.98 52.54 35.85 Further, the average performance measures over the entire dataset are summarized in Table I. The results indicate a 34% reduction in CV, 53% reduction in RMSE, and 36% improvement in motion smoothness of temporal trajectories using the present disclosure tracking framework when compared to Kinect® V1. Hence an overall improvement in the accuracy of joint motion trajectories is obtained using the present disclosure on Kinect® V1, rendering it more suitable for clinical evaluation and rehabilitation purposes.

The written description describes the subject matter herein to enable any person skilled in the art to make and use the embodiments. The scope of the subject matter embodiments is defined by the claims and may include other modifications that occur to those skilled in the art. Such other modifications are intended to be within the scope of the claims if they have similar elements that do not differ from the literal language of the claims or if they include equivalent elements with insubstantial differences from the literal language of the claims.

Various embodiments disclose methods and system for tracking body joints are able to track the trajectories of the plurality of joints associated with the plurality of body segments in an optimized manner. Here, the tracking is performed based on the direction angles and base coordinates of the plurality of initialized cylindrical models. Additionally, the optimization is performed using the state space model and the observation model. The state propagation in the state space model is determined by utilizing the state propagation rule and the weight updating is performed in the observation model by utilizing the weight updating rule.

It is, however to be understood that the scope of the protection is extended to such a program and in addition to a computer-readable means having a message therein; such computer-readable storage means contain program-code means for implementation of one or more steps of the method, when the program runs on a server or mobile device or any suitable programmable device. The hardware device can be any kind of device which can be programmed including e.g. any kind of computer like a server or a personal computer, or the like, or any combination thereof. The device may also include means which could be e.g. hardware means like e.g. an application-specific integrated circuit (ASIC), a field-programmable gate array (FPGA), or a combination of hardware and software means, e.g. an ASIC and an FPGA, or at least one microprocessor and at least one memory with software modules located therein. Thus, the means can include both hardware means and software means. The method embodiments described herein could be implemented in hardware and software. The device may also include software means. Alternatively, the embodiments may be implemented on different hardware devices, e.g. using a plurality of CPUs.

The embodiments herein can comprise hardware and software elements. The embodiments that are implemented in software include but are not limited to, firmware, resident software, microcode, etc. The functions performed by various modules described herein may be implemented in other modules or combinations of other modules. For the purposes of this description, a computer-usable or computer readable medium can be any apparatus that can comprise, store, communicate, propagate, or transport the program for use by or in connection with the instruction execution system, apparatus, or device.

The illustrated steps are set out to explain the exemplary embodiments shown, and it should be anticipated that ongoing technological development will change the manner in which particular functions are performed. These examples are presented herein for purposes of illustration, and not limitation. Further, the boundaries of the functional building blocks have been arbitrarily defined herein for the convenience of the description.

Alternative boundaries can be defined so long as the specified functions and relationships thereof are appropriately performed. Alternatives (including equivalents, extensions, variations, deviations, etc., of those described herein) will be apparent to persons skilled in the relevant art(s) based on the teachings contained herein. Such alternatives fall within the scope and spirit of the disclosed embodiments. Also, the words “comprising,” “having,” “containing,” and “including,” and other similar forms are intended to be equivalent in meaning and be open ended in that an item or items following any one of these words is not meant to be an exhaustive listing of such item or items, or meant to be limited to only the listed item or items. It must also be noted that as used herein and in the appended claims, the singular forms “a,” “an,” and “the” include plural references unless the context clearly dictates otherwise.

Furthermore, one or more computer-readable storage media may be utilized in implementing embodiments consistent with the present disclosure. A computer-readable storage medium refers to any type of physical memory on which information or data readable by a processor may be stored. Thus, a computer-readable storage medium may store instructions for execution by one or more processors, including instructions for causing the processor(s) to perform steps or stages consistent with the embodiments described herein. The term “computer-readable medium” should be understood to include tangible items and exclude carrier waves and transient signals, i.e., be non-transitory.

Examples include random access memory (RAM), read-only memory (ROM), volatile memory, nonvolatile memory, hard drives, CD ROMs, DVDs, flash drives, disks, and any other known physical storage media.

It is intended that the disclosure and examples be considered as exemplary only, with a true scope and spirit of disclosed embodiments being indicated by the following claims. 

What is claimed is:
 1. A processor implemented method for tracking body joints, the method comprising: receiving, by one or more hardware processors, motion data pertaining to a first set of motion frames from a motion sensor device, wherein the motion data comprises an initial joint motion data and a plurality of depth frames, wherein the initial joint motion data comprises locations in 3 dimensional world coordinates of a plurality of body joints associated with the plurality of depth frames and each depth frame among the plurality of depth frames comprises a set of 2D image coordinates and a depth value; initializing, by the one or more hardware processors, a plurality of cylindrical models for each motion frame among the first set of motion frames by processing the initial joint motion data and the plurality of depth frames, wherein the plurality of cylindrical models represents a plurality of body segments, connecting adjacent body joints among the plurality of body joints; and tracking, by the one or more hardware processors, the plurality of initialized cylindrical models to obtain one or more an optimized motion trajectories of the plurality of body joints in a plurality of motion frames received in succession to the first set of motion frames based on one or more direction angles and one or more base coordinates associated with each of the plurality of cylindrical models, wherein the tracking is performed by utilizing a particle filter mechanism.
 2. The processor implemented method of claim 1, wherein initializing the plurality of cylindrical models for each motion frame among the first set of motion frames comprises: obtaining a plurality of 3 dimensional point clouds corresponding to the plurality of depth frames by mapping every 2 dimensional image coordinates to the 3 dimensional world coordinates; segmenting the plurality of 3 dimensional point clouds based on a segmentation threshold to obtain a first set of segmented 3 dimensional point clouds, wherein the first set of segmented 3 dimensional point clouds corresponds to the plurality of body segments, wherein the segmentation threshold is obtained based on an Euclidean distance computed for 3 dimensional world coordinates corresponding to the plurality of body joints of the initial joint motion data; smoothening the first set of segmented 3 dimensional point clouds corresponding to the plurality of body segments prior to model fitting to reduce one or more model outliers; fitting the plurality of cylindrical models to the set of smoothened point clouds corresponding to the plurality of body segments associated with each motion frame among the first set of motion frames, wherein a set of cylinder coefficients of each cylindrical model among the plurality of cylindrical models includes radius, length, the one or more direction angles and the one or more base coordinates; initializing the radius of the plurality of cylindrical models, a mean radius value obtained from a cylinder radius corresponding to the plurality of motion frames; initializing the one or more direction angles of the plurality of cylindrical models, a mean direction value obtained from a set of cylinder directions corresponding to from the plurality of motion frames; initializing the one or more base coordinates of the plurality of cylindrical models, an average base value obtained by projecting a plurality of joint center coordinates in the initial joint motion data to an axis of a corresponding cylindrical model from the plurality of motion frames; and initializing the length of the plurality of cylindrical models, a mean length value obtained from a distance between adjacent body joints in the initial joint motion data stream for plurality of motion frames.
 3. The processor implemented method of claim 1, wherein tracking, the plurality of initialized cylindrical models to obtain an optimized motion trajectories comprises: receiving the plurality of initialized cylindrical models and a second set of segmented point cloud data corresponding to the plurality of motion frames received in succession to the first set of motion frames from the motion sensor device, wherein the initialized plurality of cylindrical models are represented as a plurality of particles, wherein each particle among the plurality of particles comprises a state represented by the one or more direction angles and the one or more base coordinates, and a weight; propagating the plurality of particles to a new state and updating the state based on a state propagation rule by utilizing a state space model, wherein the state propagation rule is based on the plurality of initialized cylindrical models; updating the weight associated with the plurality of particles based on a likelihood function and a weight update rule by utilizing an observation model, wherein the weight update rule is based on the second set of segmented point cloud data and the plurality of initialized cylindrical models; and tracking the plurality of initialized cylindrical models based on the state and the weight of the plurality of particles associated with the plurality of initialized cylindrical models to obtain the optimized motion trajectories of the plurality of body joints.
 4. A body joint tracking system, the system comprising: one or more memories comprising programmed instructions and a repository for storing the one or more motion tracking data, a 3D world coordinate and a plurality of body joints; one or more hardware processors operatively coupled to the one or more memories, wherein the one or more hardware processors are capable of executing the programmed instructions stored in the one or more memories; a motion sensor device and a joint tracking unit, wherein the joint tracking unit is configured to: receive, motion data pertaining to a first set of motion frames from a motion sensor device, wherein the motion data comprises an initial joint motion data and a plurality of depth frames, wherein the initial joint motion data comprises locations in 3 dimensional world coordinates of a plurality of body joints associated with the plurality of depth frames and each depth frame among the plurality of depth frames comprises a set of 2D image coordinates and a depth value; initialize, a plurality of cylindrical models for each motion frame among the first set of motion frames by processing the initial joint motion data and the plurality of depth frames, wherein the plurality of cylindrical models represents a plurality of body segments, connecting adjacent body joints among the plurality of body joints; and track, the plurality of initialized cylindrical models to obtain one or more optimized motion trajectories of the plurality of body joints in a plurality of motion frames received in succession to the first set of motion frames based on a set of direction angles and a set of base coordinates associated with each of the plurality of cylindrical models, wherein the tracking is performed by utilizing a particle filter mechanism.
 5. The system of claim 4, wherein the joint tracking unit is configured to initialize the plurality of cylindrical models for each motion frame among the first set of motion frames by: obtaining a plurality of 3 dimensional point clouds corresponding to the plurality of depth frames by mapping every 2 dimensional image coordinates to the 3 dimensional world coordinates; segmenting the plurality of 3 dimensional point clouds based on a segmentation threshold to obtain a first set of segmented 3 dimensional point clouds, wherein the first set of segmented 3 dimensional point clouds corresponds to the plurality of body segments, wherein the segmentation threshold is obtained based on an Euclidean distance computed for 3 dimensional world coordinates corresponding to the plurality of body joints of the initial joint motion data; smoothening the first set of segmented 3 dimensional point clouds corresponding to the plurality of body segments prior to model fitting to reduce one or more model outliers; fitting the plurality of cylindrical models to the set of smoothened point clouds corresponding to the plurality of body segments associated with each motion frame among the first set of motion frames, wherein a set of cylinder coefficients of each cylindrical model among the plurality of cylindrical models includes radius, length, one or more direction angle and one or more base coordinates; initializing the radius of the plurality of cylindrical models, a mean radius value obtained from a cylinder radius corresponding to the plurality of motion frames; initializing the one or more direction angles of the plurality of cylindrical models, a mean direction value obtained from a set of cylinder directions corresponding to from the plurality of motion frames; initializing the one or more base coordinates of the plurality of cylindrical models, an average base value obtained by projecting a plurality of joint center coordinates in the initial joint motion data to an axis of a corresponding cylindrical model from the plurality of motion frames; and initializing the length of the plurality of cylindrical models, a mean length value obtained from a distance between adjacent body joints in the initial joint motion data stream for plurality of motion frames.
 6. The system of claim 4, wherein the joint tracking unit is configured to track, the plurality of initialized cylindrical models to obtain an optimized motion trajectories by: receiving the plurality of initialized cylindrical models and a second set of segmented point cloud data corresponding to the plurality of motion frames received in succession to the first set of motion frames from the motion sensor device, wherein the initialized plurality of cylindrical models are represented as a plurality of particles, wherein each particle among the plurality of particles comprises a state represented by the one or more direction angles and the one or more base coordinates, and a weight; propagating the plurality of particles to a new state and updating the state based on a state propagation rule by utilizing a state space model, wherein the state propagation rule is based on the plurality of initialized cylindrical models; updating the weight associated with the plurality of particles based on a likelihood function and a weight update rule by utilizing an observation model, wherein the weight update rule is based on the second set of segmented point cloud data and the plurality of initialized cylindrical models; and tracking the plurality of initialized cylindrical models based on the state and the weight of the plurality of particles associated with the plurality of initialized cylindrical models to obtain the optimized motion trajectories of the plurality of body joints.
 7. One or more non-transitory machine readable information storage mediums comprising one or more instructions which when executed by one or more hardware processors causes: receiving, by one or more hardware processors, motion data pertaining to a first set of motion frames from a motion sensor device, wherein the motion data comprises an initial joint motion data and a plurality of depth frames, wherein the initial joint motion data comprises locations in 3 dimensional world coordinates of a plurality of body joints associated with the plurality of depth frames and each depth frame among the plurality of depth frames comprises a set of 2D image coordinates and a depth value; initializing, by the one or more hardware processors, a plurality of cylindrical models for each motion frame among the first set of motion frames by processing the initial joint motion data and the plurality of depth frames, wherein the plurality of cylindrical models represents a plurality of body segments, connecting adjacent body joints among the plurality of body joints; and tracking, by the one or more hardware processors, the plurality of initialized cylindrical models to obtain one or more an optimized motion trajectories of the plurality of body joints in a plurality of motion frames received in succession to the first set of motion frames based on one or more direction angles and one or more base coordinates associated with each of the plurality of cylindrical models, wherein the tracking is performed by utilizing a particle filter mechanism. 